Skewness as a probe of non-Gaussian initial conditions 
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We compute the skewness of the matter distribution arising from non-linear evolution and from 
non-Gaussian initial perturbations. We apply our result to a very generic class of models with 
non-Gaussian initial conditions and we estimate analytically the ratio between the skewness due to 
non-linear clustering and the part due to the intrinsic non-Gaussianity of the models. We finally 
extend our estimates to higher moments. 



PACS numbers: 98.80.-k, 98.65. Dx, 98.80.Cq 

The source of the initial density fluctuations which 
have led to the formation of structure, observed in the 
Universe today is unknown. Determining its nature will 
certainly be of utmost importance for the fruitful relation 
between high energy physics and cosmology. 

In models which presently attract most attention, ini- 
tial density fluctuations are generated during an infla- 
tionary phase. In the simplest inflationary models, the 
initial fluctuations obey Gaussian statistics. If this pic- 
ture is correct, the deviations from Gaussianity we ob- 
serve today were induced by nonlinear gravitational in- 
stability [0J|]- However, it is also conceivable that 
the present deviations from Gaussianity have two com- 
ponents: gravitationally induced and intrinsic, coming 
from the initial conditions rather than nonlinear dynam- 
ics . Here, we investigate to what extent an intrinsic 
component can be 'washed out' by nonlinear dynamics 
and on which scales it could be either detected or con- 
strained from above in future galaxy surveys. 

We start by deriving a general expression for the 
so-called skewness parameter, S3, including the effect 
of an initial non-Gaussianity, non-linear evolution and 
smoothing. We then estimate the normalized ./V-point 
cumulant, Sn, for a wide class of models and compare it 
with the result obtained in Gaussian models due to mild 
non-linearities . 

If the galaxies trace the spatial mass distribution, 
galaxy surveys [Q can be used to estimate the cumu- 
lants of the mass density contrast field, given by 

M N (R) = ((5 R ) N (x, Vo )) c (1) 

of the smoothed density field <5r(x, 77) = J d 3 x'Wu(\x — 
x'\)5(x', 77), where <$(x, 77) is the density field, 77 and 770 
the conformal time and its value today, and Wr is a 
window function (e.g. Gaussian or top-hat) of width R. 
The brackets in (Q) denote an ensemble average and the 
subscript c indicates that we deal with the connected 
part of the TV-point function. For a Gaussian field, all 
cumulants of order N > 2 vanish: Mn — 0. M 2 is the 
variance while M3 is a measure of the asymmetry of the 
distribution, known as skewness. We will also use the 
more common normalized cumulant, 



S N (R) = M N (R)/(M 2 (R))( N -V . 

This ratio is constant (independent of R) in the weakly 
non-linear regime |Q . To calculate the general expression 
for M 3 (i?) in the weakly nonlinear regime, we follow the 
method developed in Expanding <5(x, 77) in a pertur- 
bative series, di + 82 + 0(3) and solving the system of 
coupled Euler, Poisson and continuity equations at sec- 
ond order leads, in Fourier space, to Ai(?7, k) = Z?(?7,k) 
and 

A 2 (77, k) = (27T)" 3 / 2 / d 3 qj(q, k - q)D( V , q)D( V , k - q) 

where we consider only the fastest growing modes and 
we use the convention 

A w (77,k) = (2tt)- 3 / 2 / 6 N (r), x)e- lk x d 3 x . 

At late times where a possible source term or seed has 
decayed, the time and space dependence of the function 
D can be factorized, D(rj, k) = D + (riWk), where D + 
is the standard linear growing mode W. Perturbation 
theory gives (2) 

J(k,q) = |(1 + «) + (g/fc)Pi(/x) + §(§-«) ( 2 ) 

where the Pi is the Legendre polynomial of order I, 
fi = k • q/fcg. The quantity k is a weak function of 
fl; for fl > 0.01, k w (3/14)ft-°- 03 @. The smooth- 
ing applies order by order. In Fourier space, we have 
A_R(?7,k) = D(r],k)Wk, Wk being the Fourier transform 
of the window function. To fifth order, the skewness is 

^3 = <4,i)+ 3 (4,i^) + 0(5). (3) 

We introduce the two-, three- and four-point power spec- 
tra as (12) = -P 2 (fci)(5(ki + k 2 ), 

(123) =P 3 (k 1 ,k 2 )<5(k 1 +k 2 +k 3 ), 
(1234) c = P 4 (ki, k 2 , k 3 )5(k x + k 2 + k 3 + k 4 ). (4) 

(The Dirac 5 is a simple consequence of statistical homo- 
geneity which we assume throughout.) Here (12 . . . N) = 
(D(?7,ki)D(?7,k 2 ) . . .D(77,kjv)). The functions V 2 and 
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■p3 are also known as the power spectrum and the bispec- 
trum, respectively. Inserting the Fourier transforms of 5\ 
and 62 after smoothing in (0), expressing the correlators 
of D in terms of the power spectra (^) and performing 
one integration using the Dirac function in (^), we obtain 

M 3 (R)= I ^?V 3 (k,d)W k W q W {k+til 



(2tt) £ 



d 3 kd 3 a 

r^(*)^i(g)WjkW g W| k+q |J(k I q) + 



/ 



(2tt 
d 3 k(i 3 qd 3 p 
{2nf> 



P 4 (k, q-k, p) W 9 W p W| q+p | J(k, q-k) (5) 



For a Gaussian field, Vi — V3 — and the only non- 
vanishing contribution comes from the second term in the 
above expression. For a top hat window, this term gives 
M 3 = (34/7 -7)Mf, with 7 = -d\ogM 2 (R)/d\ogR [§. 
Note also that j(R) is the logarithmic slope of the two- 
point correlation function of the density fluctuations - 
the Fourier transform of "Piife). It is usually assumed 
that 7 > (condition of hierarchical clustering, see e.g. 

0). 

The class of models we want to analyze are those where 
fluctuations in the dark matter are induced by the energy 
and momentum of an inhomogeneously distributed com- 
ponent which contributes only a small fraction to the 
total energy momentum tensor and which interacts only 
gravitationally with the cosmic fluid. Such a component 
is denoted as 'seed' [^0|. As stressed above, we need to 
compute the A-point power spectra of the density field at 
the end of the linear regime. The comoving linear den- 
sity fluctuation D of the cosmic matter-radiation fluid 
evolves according to [pUl(| 

D + H (l - Gw + 3c 2 ) D + k 2 c 2 s D 

-§ (l + 8w ~3w 2 -6c 2 ) H 2 D = S(k,r)), (6) 

with S = (1 + w)4:TrG(fp + 3/p), f p and fp being the 
inhomogeneous energy density and pressure of the seeds. 
When the seed is a scalar field <j> with vanishing potential, 
fp+3.fp = <ft 2 - G is Newton's constant, a denotes the cos- 
mic scale factor, a dot refers to the derivative with respect 
to conformal time, H = a/a, w = P/p and c 2 = P/p are 
respectively the enthalpy and the adiabatic sound speed 
of the cosmic fluid. 

Equation (||) can be solved by a Green's function, Q, 



D(k, v ) = / g(k, v ,rf)S(k 1 ff)drf ) 
Jim 



(7) 



where rji is some early initial time deep in the radiation 
era. For the linear part of the reduced A-point function 
we then obtain 

(D(kx,ri) ■ ■ -D(k N ,rj)) c = j drji ■ ■ ■ drj N 

e(ki,7y, 771) • ••s(k iV ,?7,y7iv)(<s(ki,77i) • • • s(k N , r] N )) c . (8) 



We define the connected iV-point function of the source 
by 



(5(1) • • • S(N)) C = F N (k u ---k N ; m --- m )S (V k 



where (i) = (k j , ?7i ) . Again, the S function of the sum 
of all momenta is a consequence of the statistical homo- 
geneity. 

We now assume that the reduced Appoint function of 
the source can be replaced by its 'perfectly coherent ap- 
proximation' given by 

F N (kt, . . . ,kjv-i;v7i, . . . ,r) N ) =i sign(Fjv) x 



, kjy-i; Vi, 



, 171) • • • Fjv(ki, . . . , kjv-i; t/n, 



(9) 



(here and below, kjy is always given by k^ = — (ki + 
• • •+kjv— i))- This approximation is exact if the evolution 
equation for S is linear and the randomness is entirely 
due to initial conditions. Then the source term is of the 
form S(k,r)) = R(k)s(k, 77), where only R is a random 
variable and s is a deterministic solution to the linear 
evolution equation of S which can be taken out of the 
average (). This is the key property which renders the 
A-point function decoherent. Then can be written 
as 



Fjv(ki, . . . , kjv_i; ■ ■ ■ ,Vn) - 
S (l)--- S (A)(i?(k 1 )..-i?(k JV )) c 



(10) 



which is clearly of the form (||) . 

An important example are models with no sources but 
with non-Gaussian initial conditions for D. Such mod- 
els, like e.g. the recent x 2 Peebles model [[ll|, are always 
perfectly coherent and therefore included in our analy- 
sis: In this case D(k, rj) = R(k)d(k,rj), where R is a 
non-Gaussian random variable given by the initial con- 
dition and d is a deterministic homogeneous solution of 
Eq. Clearly, if we choose S(k,rj) — R(k)S(r/ — rj in ) 
and Q{k,rj,rf) = d(k,rf), D is of the form (0). There- 
fore, models where the non-Gaussianity is purely due 
to initial conditions are always perfectly coherent. As 
the equation of motion for D is second order, the ho- 
mogeneous solution has in principle two modes, D = 
Ri(k)di(k,?i) + R2(k)d2(k,rj), but since we shall eval- 
uate the A-point functions deeply in the matter era, the 
decaying mode will have disappeared and may thus be 
neglected in our analysis. 

Models where the source term is due to a scalar field 
which evolves linearly in time are not perfectly coherent, 
since S is given by the components of the energy momen- 
tum tensor which are quadratic in the fields. Numerical 
calculations, however, have shown that this non-linearity 
is not severe and per fect coherence is a relatively good 
approximation |l2]]l3[| . One example of this kind are ax- 
ionic seeds in pre-big bang cosmology [jl4]-[l6| for which 
decoherence has been tested and is found to be on the 
level of less than 5% for the CMB power spectrum. In 
Fig. 1 the functions D2{k,rj) and D 3 (k, k,rj) as obtained 
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by a full numerical calculation are compared to their co- 
herent approximation (j^) for the large-N limit of global 
O(N) symmetric scalar fields. This is another example 
where the scalar field evolution is linear and the only 
non-linearity in the source term is due to the e nergy m o- 
mentum tensor being quadratic in the field [pL T|J 1 3| Jl 2f| . 
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FIG. 1. The coherent approximation (dashed line) and the 
full decoherent result (solid line) for the 2- (top) and 3-point 
(bottom) functions of the large-N limit of global O(N) sym- 
metric scalar fields is shown at the end of the radiation era. 
The sign in the coherent approximation for the 3-point func- 
tion is chosen to agree with the sign for the decoherent 3-point 
function. 

For topological defects, especially for cosmic strings, 
the perfectly coherent approximation misses several im- 
portant features (like the 'smearing out' of secondary 
acoustic peaks). However, we believe that our generic 
scaling result holds also in this case, as is indicated by 
numerical simulations of global texture: even though 
global texture show considerable decoherence |jl2|| , the 
same scaling law for higher moments which we derive 
here has been discovered numerically p|. 

Under the perfectly coherent approximation Eq. (|^) 
can be factorized as the product of the N solutions, 



Dnj (ki , . . . , k^v 
term [Fjv(ki, . . 



-1,77) of the equations (0) with source 
.j-rj)] 1 ^, where kj is the 



wave number k appearing in the term c s k on the left 
hand side of (^) and the other wave numbers have to be 
considered like parameters of the source term, 



(^(k 1 , ?7 )-.. J D(k Ar ,7 7 )} c 



N 



l[D Nj {ki,...,k 



N-l 



v) 



<s(E k *) 



= Pjv(ki, ■ • ■ ,k N -i, Tjj)S(y] ki). 



(11) 



To continue, we assume that Fn is a simple power law 
in the fc, on super-Hubble scales and that it decays after 
Hubble crossing. This behavior is certainly correct for all 
examples discussed in the literature so far. We can then 
make the following ansatz 




=1 k; 



if kiT}< 1,V*G{1, 

otherwise . 



.,N} 



(12) 



Here / is a dimensionless function and kg is an arbitrary 
scale. For scale invariant seeds (e.g. topological defects) 
/ is just a constant and a — 0. For axion seeds generated 
during a pre-big bang phase, a depends on the spectral 
index of the axion field, which in turn is determined by 
the evolution law of the extra dimension Jlq| . For the 
Peebles model a is given by the power spectrum of the 
scalar field <p and / is a delta-function. Since Fjy is sym- 
metrical in the variables kj we can order them such that 
k\ > k2 > ■ ■ ■ > fcjv- 

Let us discuss the temporal behavior of the variables 
Dffj. As long as kirj < 1, the term c 2 k 2 D can be ne- 
glected in Eq. (^) and the Green's function is a power 
law. At fci?7 ~ 1 the source term decays and as long as 
a perturbation remains super horizon, it just grows like 
rj 2 , so that for kjf] < 1 < k\r), 



D Nj w 5(l//ci)fc 1 



2+3/N , 



where 



1(V) 



4ttG 



V 



2-3/N 



" G{^i)f(i)^-^ d 4 , 

in V 



and we have to take the part of the integral above which 
converges when rji n — > 0. 

Once the perturbation enters the horizon it either 
starts oscillating with roughly constant amplitude or con- 
tinues to grow oc rf , depending on whether kj enters dur- 
ing the radiation or matter dominated era. At late time, 
V ^ Veq aud kr] ^> 1, we therefore obtain 



D Nj « 5 (l/fc 1 )fc- 2+3/Ar (A :i /fc J ) 

N 



JL. 

>Je„ J ^ "I ' le 1 

2 



(rjkj 



if kjT) eq > 1 

if kjT] eq < 1 
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where rj eq is the time of equality between the matter and 
radiation densities. Defining < j eq < N so that kjrj eq > 
1 for all j < j eq we obtain for the connected ./V-point 
function 



Pjv(ki,. 



V/V-l 



n 



n 



kji] e 



(13) 



Using this result for the ordinary power spectrum, V2, 
we can express Vn is terms of products of V% as 



•pjv(ki,---,kjv_i,7?) 

N 

■3(l-N/2) 



3/2 ' 



3/2 



(14) 



For the class of models considered and under the assump- 
tion of perfect coherence, we have determined the con- 
nected A-point power spectra in the linear regime which 
are the input of the skewness (||). 

M3 has two contributions: A linear one due to the 
initial non-Gaussianity (contained in V3) and one due 
to non-linear clustering which induces skewness even in 
an originally Gaussian distribution of perturbations; it 
contains a Gaussian part ("Pf ) and a non-Gaussian term 
(Pi)- We decompose the skewness as 

We want to estimate the ratio of these two contributions. 
Under our approximation (fbf), the first term of re- 
duces to 



d 3 kd 3 k' 



W k W q W lu+k ,^p2(k)V 2 {q)p2{\^ + q\) 



1,-3/2 



g(l/fc m ax) 3 (fcg|k + q|) 3 / 2 
g(l/k)g(l/q)g(l/\k + c l \)k 9 Jl 



(15) 



where fc max = max{/c, q, |k + q|}. M 3 is given by the 
second and third terms in (^|). 

To estimate analytically the ratio /M^ NL ^ 



si L) /s 3 



(NL) 



, we make the following approximations: 



We assume that V2 is a simple power law within the 
range of scales of interest, namely all the modes 
which enter the horizon during the radiation era, 
this is O.l/i^Mpc^ 27r/fc < 20ft, _2 Mpc, name. 
V 2 (k) = k- 3 (k/K)'' . 

We also assume that g(rj) oc rf . 

We replace the window function by a simple cut-off 
at k = l/R. 

For symmetry reasons we may integrate over the 
triangle q < k < R and then multiply the result by 
2. 



- Since in our integration region, q < k, we replace 
|k + q| by k. 

With these approximations the angular dependence of 
the integrand disappears and the integrals over k and q 
in (^|) can be trivially performed leading to 



Mi L \R) ~ 



M. 



(NL) 



for 



4(k,R)-^/ 2 
(27r) 4 3 7 (3 + 7/2 + r) 
7 > and 3 + 7/2 + r > 

(KR)-^ r 

7 > 0, 



(16) 



(2tt) 4 7 2 



for 



where we have just considered the Gaussian contribution, 
P 2 2 toMf L) . 

Since fc» is just the scale beyond which the density 
contrast (D(x) 2 ) R= i/ k ~ V 2 {k)k 3 is larger than unity 
and non-linearities become important, we define the non- 
linearity scale i?ii n = 1/fc* ■ The ratio between the skew- 
ness due to the non-Gaussianity in the linear perturba- 
tion and the one due to dynamical nonlinearities is then 



5'. 



(L) 



47 



Si NL ) 3(3 + 7/2 -r) 



R 



7/2 



(17) 



This is our main result. It is readily checked that the 
non-Gaussian contribution, Va, to M^ NL ^ behaves just 
like the contribution M^ NL ^ and thus only modifies the 
pre-factor in , which should not be taken too seriously 
in view of the relatively crude approximations which we 
have employed to obtain our result. 

This computation of the skewness is easily generalized 
to higher moments. As our computation shows, linear 
non-Gaussianities scale like 



M { N L) (R) oc {R/R^)- N ^ 2 



(18) 



The dominant non-linear contribution to the connected 
iV-point function which is also present in Gaussian theo- 
ries contains N — 2 second order terms D 2 and there- 
fore scales like 



M 



(NL, Gauss) 
N 



(R) cx (R/R Un )-( N -^ 



(19) 



The lowest order non-linearity for a generic non-Gaussian 
model, however just comes from the non-Gaussian term 
with N + 1 factors of D. The non-Gaussian non-linear 
corrections therefore generically scale like 



M 



(7VL,noGauss) 
N 



(R) oc (i?/ J R lin )-( JV + 1 )^ 2 



(20) 



Only for N — 3 the two terms ( ]19| ) and ([20j) scale in the 
same way. For all higher A's the non-Gaussian contribu- 
tion dominates in the mildly non-linear regime, R > i?n n . 
From Eq. (|2(]) we infer that in on large scales the ratios 
for all reduced A-point functions very generically scale 
like 



4 



S ( N L) (R) 

s { N NL \R) 



oc 



R 

Rlin 



7/2 



(21) 



This expression agrees with other analytic predictions 
as well as numerical simulations in a global texture 
model ||. The agreement with the texture simulations 
which are decoherent suggests that the validity of our re- 
sult extends beyond the conditions under which Eq. j2l| ) 
was derived. More important than decoherence is that 
the source term decays at late times and therefore the 
density perturbations just evolve according to the ho- 
mogeneous solution. This implies that at late times the 
iV-point functions behave like the homogeneous growing 
mode to the iVth power, while the reduced TV-point func- 
tion induced by non-linear clustering from Gaussian per- 
turbations scales like the growing mode to the 2( J ?V — l)th 
power. Since topological defect sources decay on sub- 
horizon scales, we conclude that the derived scaling be- 
havior is also valid for them (this argument will be ex- 
panded in our follow up publication ||). 

Our result implies that on small scales (R < Ru n ), the 
dominant contribution to the cumulants comes from non- 
linear Newtonian gravitational clustering, and the Gaus- 
sian term actually dominates. Intrinsic deviations from 
Gaussianity are difficult to detect on small scales. Hence, 
we should look for signs of intrinsic non-Gaussianity at 
large scales (R > i?i; n ). This suggestion was expressed 
earlier based on qualitative physical arguments j|; how- 
ever, our present result is derived from first principles for 
a specific class of initial conditions - coherent seeds. 

If galaxies trace mass, the measurements of the two- 
point correlation function suggest Ru n ~ 10/i _1 Mpc and 
j(R) « 1.8 for lOkpc < hR < 15 Mpc (here h is the usual 
parameterization for the Hubble constant in units of 100 
km s _1 Mpc _1 ); the slope 7 becomes steeper at larger 
separations R [[T|Jlg|] . A frequently considered theoretical 
possibility for long- wave tail of the initial V2 (k), called 
the Zel'dovich-Harrison spectrum, would give 7 = 4 at 
large separations. Hence, we can expect all SnS to "blow 
up" with increasing scale for the class of non-Gaussian 
models considered here, in contrast with models with 
Gaussian initial conditions. The available measurements 
of Ss(R) and Si(R) do not show such a rise with scale 
and have already been used to constrain texture mod- 
els ||. Likewise, there are indications that the existent 
data from the APM Galaxy Survey may may already ex- 
tend to sufficiently large scales to constrain the \ 2 Pee- 
bles model [|lQf| . With surveys presently underway like 
the Sloan Digital Sky Survey |^(|, the prospects for us- 
ing the approach outlined here to probe the statistics of 
the cosmological initial conditions will become even bet- 
ter. 

In this work we derived a scaling law for the "intrinsic 
to induced" skewness ratio ( |i~7|) for coherent seeds. We 
also showed how to generalize this law to higher cumu- 
lants. We plan to follow these calculations with more 
detailed predictions for coherent seed models and to con- 
front our analytic results with numerical simulations as 



well as observational data from galaxy surveys j|. Let 
us also repeat that the derived scaling laws seem to be 
more general than their derivation as they have been 
obtained numerically for global texture which are deco- 
herent seeds. We actually believe that the origin of the 
scaling laws is not coherence but mainly the decay of the 
sources at late time and we therefore conjecture that they 
hold also for topological defects. 
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